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Probing the physical and chemical structure of the CS core in 
LDN 673. Multitransitional and continuum observations 
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ABSTRACT 

High-angular resolution observations of dense molecular cores show that these cores can be 
clumpier at smaller scales, and that some of these clumps can also be unbound or transient. 
The use of chemical models of the evolution of the molecular gas provides a way to probe 
the physical properties of the clouds. We study the properties of the clump and inter-clump 
medium in the starless CS core in LDN 673 by carrying out a molecular line survey with 
the IRAM 30-m telescope toward two clumps and two inter-clump positions. We also ob- 
served the 1.2-mm continuum with the MAMBO-II bolometer at IRAM. The dust contin- 
uum map shows four condensations, three of them centrally peaked, coinciding with pre- 
viously identified sub-millimetre sources. We confirm that the denser clump of the region, 
n ~ 3.6 x 10 5 cm -3 , is also the more chemically evolved, and it could still undergo fur- 
ther fragmentation. The inter-clump medium positions are denser than previously expected, 
likely n ~ 1 x 10 3 — 1 xlO 4 cm -3 due to contamination, and are chemically young, similar 
to the gas in the lower density clump position. We argue that the density contrast between 
these positions and their general young chemical age would support the existence of transient 
clumps in the lower density material of the core. We were also able to find reasonable fits 
of the observationally derived chemical abundances to models of the chemistry of transient 
clumps. 
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1 INTRODUCTION 

It has been long known tha t molecular clouds are highly structured 
(e. g.. lBlitz& Starklll986h and their structure is grea tly affected 
by the motions induced by supersonic turbulence (e. g.. lScalo et al.l 
1 19981). self-gravity of the gas and magnetic fields i nside the clouds 
dMcKee & Ostriked[2007l ; iKainulainen et alj|200gj) . All these pro- 
cesses control the formation and evolution of the density enhance- 
ments, of different scale sizes and densities, such as the cores and 
clumps that will finally give birth to stars. But there are still great 
uncertainti es to identify the connection between sta rless cores and 
protostars djohnstone et al j200CtlSmithetal.l2008» . Starless cores 
are not all the same, despite their overall similarity in structure 
dKeto & Caselli[2~008h . and differences in total mass, density and 
temperature might account for the differences in dyna mical proper- 
ties, structure, and future ev olution of starless cores dKeto & Field! 
l2005tlKeto & Caselll2008T) . 

Higher-angular resolution observations are also finding that 
dense cores in molecular clouds that appeared homogeneous in 



single-dish observations, are clumpier at smaller scales dPeng et al .1 
ll998l : lMorata et al]|2003l) , showing structures as small as 0.02 pc, 
and masses as low as O.OIMq. Many of these smaller clumps are 
unbound and/or showing evidence of being transient dPeng et al.l 
ll998l : lMorata et alj2005t) . and will never be able to form low-mass 
stars or even brown dwarfs. The mix of bound and unbound struc- 
tures in dense cores is also found in several regions, such as the 
Pipe Nebula, where recent molecular line and continuum observa- 
tions (iLombardi et al. 2006; iMuench et al.l 2007: Rat hborne et all 
l2008t iFrau et al .1120100 found numerous cores more than 100, most 
of whi ch appear to be p ressure confined, and gravitationally un- 
bound dLada et all2008h . 

Observations of the emission of molecular lines at millimetre 
and sub-millimetre wavelengths in cloud cores combined with the 
modelling of the chemistry of the gas provides a way of obtaining 
information on the physical structure and the chemical and physi- 
cal evolutionary stages of the cores. We proposed a time-dependent 
chemical model that also explored the consequences of the pres- 
ence of unresolved and transient structures in the gas that would 
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form and disperse in a timescale of ~ 1-2 Myrs dTavlor et al.l 
1 19961) , in order to explain the systematic differences between CS 
and NH3 lines dPastor et al.ll99ll:lMorata et aljl997h. Simulations 
of the evolution of these cores 1 Garrod et al.ll2005 , 20061) find that 
clouds that are ensembles of such transients have a clearly differ- 
ent chemistry from a 'traditional' static cloud. The gas chemistry 
appears to be "young" at all times, and the re-cycling of the mate- 
rial frozen out onto dust grains produces a general molecular en- 
richment of the clouds, even after re-expansion of the transient 
structures. The background gas in which these inhomogeneities 
are embedded would be fairly diffuse, but chemically enriched. 
These chemical enhancements might also account for the variety 
of chemistries observed in diffuse clouds. 

We carried out interferometric high-angular res olution obser- 
vation s of several molecules (CS, HCO + and N2H + ) dMorata et al .1 
2003, hereafter MGE03), which we later combined with singl e-dish 
intermediate-angular resolution maps dMorata et al.l I2005L from 
now on MGE05) towards the starless CS core in LDN 673 (d = 
300 pc), in order to test the predictions of the chemical models. 
The combined single-dish and interferometer maps showed emis- 
sion of both background and clumped gas, with a clear segrega- 
tion of clump properties between the northern and southern halves 
of our observed region, and allowed us to identify 15 resolved 
clumps in our data cube. The derived clump masses are well be- 
low the virial mass, which would point to their being transient, ex- 
cept for the more massive one, which might have a mass, ~ IMq, 
closer to the virial mass. The starless core appears to be consti- 
tuted by a heterogeneous medium of condensations, of various 
densities and at different stages of chemical evolution, in agree- 
ment with theoretical studies that postulate the existence of tran- 
sient clumps or the transient nature of dense core s generated by dy- 



Vazquez-Semadeni et al] 120051 : IVan Loo et al.l l2008h . Re 



WhvattetalJfeOlCh also found evidence of a heterogeneous 



namical flows within molecular clouds ( seee. g., Falle & Hartquistl 
120021 ; 
cently, 

medium in scales of less than 0. 1 pc near HH objects, as traced by 
strong HCO+ (3-2) emission. These clumps would have gas vol- 
ume densities > 3 x 10 4 cm" 3 . 

In order to study the properties of the clump and inter-clump 
gas in the starless CS core in LDN 673, we selected two positions 
associated with identified clumps (CL1 and CL6) and two posi- 
tions where the inter-clump gas would be dominant (where we did 
not detect any clump). A multitransitional survey of several early- 
and late- type molecules in these positions allows us to sample the 
chemical composition of the gas and compare it to the predicted dif- 
ferent chemistry of the pre- and post-clump gas in the models. We 
additionally observed the dust continuum emission in LDN 673. 
The structure of this paper is as follows: in Sect. [2] we describe the 
IRAM 30-m spectral line and continuum observations. In Sect. [3] 
we describe the characteristics of the detected spectra and of the 
dust continuum emission. The analysis of the observational results 
and the determination of the physical parameters of the gas and dust 
are shown in Sect. [4] Finally, Sect.[5]contains the discussion of the 
results of our analysis and how they can be related to the previous 
observations, the chemistry of the clouds and the structure of the 
core. 



2 OBSERVATIONS 

The spectral line observations were carried out in 2005 August us- 
ing the 30-m IRAM telescope in Granada (Spain). We used the ca- 
pability of the ABCD multi-receiver system to observe 10 different 



Table 1. Lines observed with the 30-m IRAM telescope 



Transition 



CCH 1-0 |-| F=2-l 
CCH 2-1 |-| F=3-2 
CCH 3-2 |-| F=4-3 



CN 1-0 



3 5 



CN 2-1 ^ i 



1 3 

2 • 2 



2 ' 2 2 ' 2 
CS 5-4 

C-C3H2 2i j2 -lo,l 
C-C3H2 4 M -3 ,3 
HCN 3-2 

H 2 C0 2i, 2 -li,i 
H 2 CO 3i, 3 -2i,2 
H 2 CO 3i,2-2i,i 

H 13 CO+ 1-0 
H 13 CO+ 2-1 
H 1:i CO+ 3-2 

no §,§-§.§ n+ 
no - --- - n+ 

2' 2 2*2 

SO 3 2 -2i 
SO 4 5 -3 4 
SO 65-54 
SO 7 6 -6 5 

S0 2 31,3-20,2 
SO2 7l,7-6o,6 



Frequency 


HPBW 


-Boff 


(GHz) 


(arcsec) 




87.316925 


28 


0.78 


174.663222 


14 


0.64 


262.004260 


9 


0.46 


113.490982 


22 


0.74 


226.874764 


11 


0.54 


244.935606 


10 


0.50 


85.338906 


29 


0.78 


150.851899 


16 


0.68 


265 886432 


9 


0.45 


140.839515 


18 


0.70 


111 1 1 1 A A Q 
Zl L.ZL 1445 


1 1 


0.57 


iiJ.U? I 1 I J) 


I ] 


0.54 


86.754330 


28 


0.78 


173.506782 


14 


0.64 


260.255480 


10 


0.46 


150.176459 


16 


0.68 


250.436845 


10 


0.48 


99.299905 


25 


0.76 


206.176062 


12 


0.58 


219.949433 


11 


0.55 


261.843756 


10 


0.46 


104.029410 


24 


0.76 


165.225436 


15 


0.66 



Table 2. Coordinates of the positions selected for the spectral line observa- 
tion. 



Position 


Counterpart 


R.A.(J2000) 


Dec. (J2000) 


CL1 


N 2 H+ peak 


19:20:51.747 


11:13:49.50 


CL6 


CS peak 


19:20:50.003 


11:14:53.00 


ICLN 


Inter-clump N 


19:20:51.701 


11:15:30.00 


ICLS 


Inter-clump S 


19:20:54.501 


11:14:10.00 



molecules (and a total of 23 transitions) with just 6 frequency se- 
tups in the 3, 2, and 1-mm bands. TableQ]shows the transitions and 
frequencies observed. We used the VESPA autocorrelator as a spec- 
tral backend, which provided a total bandwidth of 80 MHz, and se- 
lected a 20 kHz channel spacing for the receiver at 100 GHz, and a 
40 kHz channel spacing for the other three receivers. The achieved 
velocity resolutions range from 0.4 to 0.07 km s _1 from the 1 to the 
3-mm bands. The main-beam efficiencies and the half-power beam 
widths at the observed frequencies are also listed in Table Q] We 
used the frequency-switching mode with a 7.9 MHz throw, except 
for the configuration that observed simultaneously the CCH (1-0), 
C-C3H2 (4i,4-3 , 3 ), H 2 CO (3i, 2 -2i,i), and SO (7 6 -6 5 ) transitions, 
where we used a 15.8 MHz throw. We obtained system tempera- 
tures, in T* scale, of 85-160 K at 100 GHz, 215^165 K at 150 GHz, 
230-320 K at 230 GHz, and 470-680 K at 270 GHz. We used the 
GILDAsQ package of IRAM to reduce, analyse, and display the 
spectral data. 



1 |http: //www.iram. f r/IRAMFR/GILDAS| 
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19*20 m 55 - 50" 19*2<f55' 50' 

RA (J2000) 



Figure 1. Positions observed with the 30-m telescope in the spectral line 
observation (crosses) overlapped with the greyscale and contour channel 
map of the CS (2-1) emission obtained by MGE05. The two concentric 
circles mark the approximate 30-m beam sizes at 1 and 3 mm, respectively. 

Table [2] gives the coordinates and names of the four observed 
positions. Figure[T]plots the location of the observed positions over 
the combined BIMA-FCRAO map of the CS (J=2-5-l) emission of 
MGE05. The four positions were selected according to the follow- 
ing criteria: (1) two positions centred in two of the CS clumps de- 
tected by MGE05: CL1, the most massive and chemically evolved 
clump, and CL6, a chemically young clump found in the north- 
ern region of the BIMA map; (2) two inter-clump positions, one 
in the northern region (ICLN) and another in the southern region 
(ICLS), selected from the study of dumpiness of the CS emission 
in MGE05. 

The continuum observations were carried out in February 
2007 with the 30-m telescope using the MAMBO-II 117-channel 
bolometer at 1.2 mm. The bolometer effective frequency is 250 
GHz and the half-power beam-width at 1.2-mm is 11". The source 
was observed using the On-the-fly mode, with a secondary chop- 
ping of 42". The telescope was scanning in azimuth at a speed 
of 8" s . Total observation time was 1.4 h, with 57 minutes on 
source. During the observation, the zenith atmospheric opacity at 
225 GHZ was 0.32. The resulting rms was ~ 6 mjy beam -1 . 
The map was centred at a(J2000) = 19 h 20 m 52!20; <5(J2000) = 
11° 14'45"5 and it covers an effective area of about 9!0 x 6^5 along 
the NE-SW direction. We used the MOPSI C package of IRA M to 
reduce the continuum data and the MIRIAD iSault et aljl995h and 
GILDAS software package to analyse and display the continuum 
data. 



3 RESULTS 

3.1 Spectral line observations 

Figures [2] and [3] show the spectra obtained in the four selected po- 
sitions for the detected molecules (19 of the 23 lines observed). 
Tables lATI to I A4I show the line parameters obtained using a Gaus- 
sian fit for the transitions detected in the CL1, CL6, ICLN, and 
ICLS positions, respectively. Table |A5l shows the hyperfine struc- 
ture fit parameters (obtained using the HFS method of the CLASS 
package) for the detected CN, CCH and NO transitions (Figure [4] 




468 10 468 10 468 10 468 10 
v lsr (km/s) v^b (km/s) v LSR (km/s) v LSR (km/s) 



Figure 2. Spectra of the detected transitions of the H2CO, H w CO+, CN, 
C-C3H2, and CCH molecules at the four selected positions in the CS core 
of LDN 673, in the velocity range 3-11 km s _1 . The dotted vertical lines 
indicate the Vlsr °f me H2CO (3i 3— 2i 2) transition for the respective 
position. 



shows the full spectra of all the detected hyperfine structure lines). 
Finally, Table |A6l gives the upper limits for the line intensities of 
the transitions not detected at each position. The CCH (3—5-2), CS 
(5 — 5-4), SO (76—5-65), and SO2 7i,7-6o,6 lines were not detected at 
any position. The first three lie in the 1-mm band, which presents 
higher T sys , and typically shows a higher rms than the 3- and 2-mm 
band values in our spectra. 

All the 19 detected transitions are at the CL1 position, whereas 
the other three positions show detection in 10 (CL6 and ICLN) or 
1 1 (ICLS) transitions. If we compare the intensity for the same line 
across the four positions, the emission is always more intense at 
the CL1 position. Lines at the CL6 position tend to be more intense 
than at ICLS, while lines at the ICLN position are typically the less 
intense ones. 

There are not important differences in the line central ve- 
locities among the four positions. The four observed positions 
show i>lsr velocity differences of 0.1-0.2 km s" 1 between 
them, and are in agreement with the velocity pattern found by 
MGE05. Linewidths toward CL1 tend to be the narrowest (~ 0.4- 
0.5 km s~ , typically). This suggests that the turbulence is lower in 
CL1 than in the rest of the region. 
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Figure 3. Same as Fig.[2]for the SO, S0 2 , NO, and HCN molecules. 



3.2 Dust continuum emission 

Figure [5] shows the 1.2-mm continuum map obtained with the 
MAMBO-II bolometer. We find four emission condensations with 
different degrees of extension and intensity. T hese four condensa - 
tions were previously detected at 850 by IVisser et"al I J2002h . 
Table [4] gives the emission parameters derived for the se four con- 
densat ions and the association with the sources of IVisser et"al] 
d2002h. The s e four condensations were also recently identified by 
iTsitali et al.1 ( bold) using the Spitzer Space Telescope, which they 
classify as starless cores. 

The three southern condensations show an approximately 
rounded shape and with centrally peaked emission, whereas the 
northern condensation rather looks as an extended patch of weak 
emission. The most intense emission peak (the central clump in 
Fig.|5j. SMM5, coincides with CLl and with the N2H 4 " emission of 
MGE03. Another, slightly smaller, bright condensation is SMM3, 
with a similar value of the emission peak and located about ~ 50" 
west of SMM5. A third condensation, again slightly smaller, but 
with a similar emission peak to the previous two, is SMM8, located 
~ 2.' 5 SE of SMM5, near the edge of the 1.2 mm observation's 
field of view. This condensation was also detected in H 13 CO + 
(1-0) and C 34 S (2-1) in single-dish observations (MGE05). The 
fourth condensation, associated with SMM6, presents diffuse and 
extended dust emission in an approximately E-W strip ~ 2' north 
of SMM5 that coincides with a ridge of molecular emission found 
in an E-W direction in the maps obtained with the FCRAO tele- 
scope (MGE05). 

Among the four positions observed in spectral lines, only CLl 
is associated with a dust peak matching very well with SMM5. CL6 
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Figure 4. Spectra of the detected lines showing hyperfine structure: a) CN, 
CCH and NO lines at 3-mm; b) CCH (2-1) and CN (2-1) lines in CLl. 
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LDN 673 
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Figure 5. 1.2-mm dust continuum map obtained with the MAMBO-II 
bolometer after convolving with a 15" Gaussian beam. The nns noise of 
the map is 7 mJy beam -1 . The minimum contour is 21 mjy beam -1 , with 
increments of 7 mjy beam - 1 . The circles show the beam size of the 30-m 
telescope at 3mm, ~ 25", centred at each of the positions of the spectral 
line survey (solid squa res). The solid tria ngles mark the positions of the 
SMM sources found bv lVisser et aI1fe002l) . 
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and ICLN are associated with weak dust emission at 2<r, and ICLS 
is associated with even weaker emission. 



4 ANALYSIS 

4.1 Radiative transfer analysis 

The detection of multiple transitions at the same position of some 
molecules allows us to constrain the physical properties of the 
emitting gas. We used th e RADEX radiative transfer model code 
jVan der Tak et al.ll2007l) . which treats optical depth effects with 
an escape probability method, to find the set of molecular column 
densities, N(X), volume densities, n(H2), and kinetical temper- 
atures, Tk, that best reproduced the observed lines. Appendix |B1 
gives more details on how the calculations were done. 



f 5 
6 



0,= 5O", N(H 2 C0)=2.Bxl0 cm 



e„=1000". J»(H 2 C0)=2.5xl0 13 cm ; 




= 1000". J»(H 2 C0)=3.0xl0 13 cm 1 



10 15 20 25 10 15 20 25 30 

T k (K) 



4.1.1 Results of the RADEX analysis: uncertainties 

Table [4] shows the best-fit results we obtained from our xt anal- 
ysis using RADEX on the set of lines from the H2CO, SO, and 
H 13 CO + molecules, for which at least 3 transitions have been de- 
tected at each position. Table[5]shows the best set of parameters that 
we found for C-C3H2, HCN, and SO2, assuming the temperature, 
volume density and size derived from the previous best fit, plus the 
results from the hyperfine-line fitting for CCH, CN, and NO. 

We estimate that the volume density determination has an un- 
certainty of 50-70% of the best fit value (shown in Table 4) at the 1- 
sigma level. The column densities, N(X), are less well constrained 
than the volume density. The derived column densities of all the 
molecules have an uncertainty of ~ 65% on average from the best- 
fit value for the CL1 position. For the other positions, the column 
densities are only constrained within a factor of 2 of the best value. 

We alternatively explored how well constrained were the re- 
sults of our analysis by fixing the H2CO column density and the 
filling-factor for each position to the best fit values, and running a 
range of models with different kinetical temperatures and volume 
densities. Figure [6] shows the goodness-of-fit over the parameter 
space spanned by n(Ha) and Tk. Similarly to the results in Ta- 
ble [4] n(H2) is relatively well-constrained for all positions, but, 
except for CL1, that is not the case for Tk. Temperatures are rather 
badly determined for the two inter-clump positions, probably due 
to the lower volume density and the lower SNR of the 1-mm lines. 
The best fit values of Tk show clear differences between the CL1 
and CL6 positions, Tk — 10 K, and the inter-clump positions, 
Tk — 23 K (see also Sect. 14. 1 ."2b . However, the uncertainties in 
the inter-clump positions are too large to be significant. We adopt a 
value of 20 K for the inter-clump gas, since their lower visual ex- 
tinct ion may allow the interstell ar UV radiation to heat the gas (see 
e.g jNeiad & Wagenblastll 19991) . 

As shown in Fig. [TJ the 2 and 3 mm beams of the telescope 
at both inter-clump positions could have contribution from higher 
density gas. This effect could be important for the ICLS position, 
which is the one we studied more carefully (the effect of the emis- 
sion of CL6 on ICLN is probably very small because the emission 
in CL6 is much less intense). We used the 1.2-mm continuum map, 
with an angular resolution ~ 10'.' 5, to estimate the contribution to 
the intensity of the lines at the ICLS position coming from emis- 
sion from the CL1 position. In order to have an upper limit estimate 
of the error beam of the telescope at 2 and 3 mm, we convolved 
the 1.2-mm continuum map with a Gaussian with the same HPBW 
as the main beam of the telescope at 2 and 3 mm obtained from 



Figure 6. Contours of constant xt value (the goodness-of-fit) over the pa- 
rameter space spanned by the volume density, n(H2), and Tk for fixed 
values of the H2CO column density, N (H2CO), for each of the observed 
positions. The contours mark the confidence levels of 68.3%, 90.0%, and 
99.9%. Red (grey) crosses mark the best-fit parameters for each position. 

iGreve et al.l i ll 9981) and lBensch et all ( l200lT) . Then, we compared 
the intensities at the clump and inter-clump positions for the origi- 
nal and the convolved maps. We calculated the contribution of the 
error beam at the CL1 position to the H 2 CO, SO, and H 13 CO + 
lines at the ICLS position as representatives of the different beams 
we used. We found that the error beam contribution to the molecu- 
lar line emission in ICLS can be < 0.1-0.5 K, or < 18-50% range 
of the intensity of the detected lines. Thus, though not negligible, 
we expect a minor contribution of the higher density gas to the lines 
in the ICLS position. In order to check how much the results of 
the analysis for the ICLS position can be affected by this contribu- 
tion, we remade the analysis after subtracting the estimated inten- 
sity coming from the CL1 position in the lines measured in ICLS. 
The resulting number density for ICLS becomes ~ 3 times smaller, 
~ 6.6 x 10 3 cm~ 3 , while the column densities of H 2 CO and SO 
are ~ 2 times larger, and H 13 CO + remains practically constant. 

4.7.2 Results of the RADEX analysis: volume density, excitation 
temperature and line opacity 

The line intensities and line ratios of our data are best reproduced if 
we use beam-filling factor and kinetical temperature, r\\ j < 1 and 
T K ~ 10 K, for CL1 and CL6, but r& f ~ 1 and T K ~ 20 K give 
better results for the two "inter-clump" positions. We performed 
an additional test to check if the "best-values" of the filling-factors 
are consistent with previous data. From the maps by MGE05, we 
calculated the main beam antenna temperatures at each of the four 
positions after convolving the maps with beams of FWHM 14" and 
28" and tried to estimate the size of the "apparent" emitting region. 
The results for the ICLN and ICLS positions agree with r\ b j ~ 
1. For CL1 and CL6, the results agree with S ~ 30" and 45", 
respectively, not far from the adopted value 6 a ~ 50". 

The volume densities yielded by the \t analysis show a clear 
distinction between the physical characteristics of the four posi- 
tions traced by the molecular line emission. The volume density is 
highest in CL1, ~ 3.6 x 10 5 cm~ 3 , ~ 5 times larger than the den- 
sity derived for CL6, and more than one order of magnitude larger 
than for the inter-clump positions. 
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3. Parameters of the condensations found in the 1.2-mm dust emission 



Peak Flux" 





Position 


intensity 


density 


FWHM 6 


Mass c 


N(H 2 ) d 


n(H 2 ) d 




No. 


R.A.(J2000) 


DEC.(J200) 


( Jy beam ) 


(Jy) 


(arcsec) 


(M ) 


(10 22 cm- 2 ) 


(10 5 cm~ 3 ) 


Association e 


1 


19:20:46.0 


11:16:20.0 


0.038 


0.436 


75 


3.9 


1.6 


0.7 


SMM6 


2 


19:20:48.2 


11:14:10.5 


0.066 


0.253 


41 


2.9 


4.0 


3.3 


SMM3 


3 


19:20:52.0 


11:13:49.5 


0.066 


0.339 


49 


5.1 


4.9 


3.4 


CL1/SMM5 


4 


19:21:02.4 


11:11:43.5 


0.065 


0.186 


35 


1.6 


3.0 


2.8 


SMM8 



a Flux density integrated above the 2<r level after convolving with a 15" Gaussian beam. 

b FWHM of the dust emission condensation, %whm = sjQ's ~ ®b' wnere @b ls me beam-size and 9 S = 2^J A/it. A is the area on the map of 
the pixels with intensity over half of the local continuum emission peak. 

c Total mass calculated from the flux density integrated above the 2<r level, using the expression from lFrau et al.|j2010h . for a distance of 300 pc. 
d iV(H2) andn(Fi2) are the FWHM averaged values of column and volume densities, respectively, calculated using the expressions from 
iFrau et al assuming dust temperature, T(j uat = 10 K, a standard dust-to-gas ratio of 100, a dust absorption coefficient at a frequency 

of 250 GHz, K25Q = 0.0065 cm 2 g —1 (tak en as a medium value between dust grains with thin and thick ice mantles for a volume density of 
~ 10° cm" - 3 . lQssenkopf & Hennin j 119941) , which should be typical for starless cores. 
e SMM: sub-millimetre sources from lvisser et al . (2002). 

Table 4. Physical parameters derived from the best-fit model defined by the minimum value of the xi function of a series of RADEX models compared to the 
lines of H2CO, SO, and H 13 CO+ observed in the four positions. 



Position 


Species 


6 S 
<") 


Tk 
(K) 


n(H 2 ) 
(10 5 cm -3 ) 


N(X) a 
(cm- 2 ) 


Transition 


r r 1 b 
± ex 

(K) 




CL1 


H 2 CO 


50 


10 


3.6 


8.4 X 10 12 


2i,2-ll,l 


6.3 


1.0 




SO 








4.8 x 10 13 


3 2 -2i 


8.2 


2.0 




H 1:i CO+ 








8.3 x 10 11 


1-0 


9.5 


0.4 


CL6 


H 2 CO 


50 


10 


0.7 


2.8 x 10 13 


2l,2-ll,l 


4.7 


4.4 




SO 








1.3 x 10 14 


3 2 -2i 


6.7 


7.2 




H 13 CO+ 








4.6 x 10 11 


1-0 


5.0 


0.4 


ICLN 


H 2 CO 


1000 


20 


0.2 


3.0 x 10 13 


2i,2-li,l 


3.9 


3.6 




SO 








6.6 x 10 13 


3 2 -2i 


4.9 


2.9 




H 1:i CO+ 








4.0 x 10 11 


1-0 


4.2 


0.3 


ICLS 


H 2 CO 


1000 


20 


0.2 


2.5 x 10 13 


2i,2-ll,l 


4.0 


6.1 




SO 








7.6 x 10 13 


3 2 -2i 


5.0 


4.9 




H 13 CO+ 








4.0 x 10 11 


1-0 


3.8 


0.5 



a Beam-averaged column densities. 

b Calculated for the lower frequency transition (the best determined in the observations). 



Table[4]also shows the values of excitation temperatures, T ex , 
and line opacities, r, yielded by the "best-fit" models for each po- 
sition and molecule: H 2 CO (2 1)2 -1 M ), SO (3 2 -2i), and H 13 CO+ 
(1-0). We can see that the H2CO and SO lines tend to be optically 
thick or very thick, especially for the CL6, ICLN, and ICLS posi- 
tions, which is also reflected in apparent larger values of the column 
density, compared to the ones found at CL1. This effect is possibly 
due to the lower volume density in CL6, ICLN, and ICLS. The ro- 
tational levels involved in the H2CO and SO lines are not of the 
fundamental rotational level, which makes them more difficult to 
populate at lower densities. In order to reproduce the observed line 
intensities, we need to compensate by using higher column den- 
sities (and line opacities). This effect is not seen for the H 13 CO + 
lines, involving the fundamental rotational level, which are rela- 
tively optically thin and with CL1 column densities larger or simi- 
lar to the ones found at the other positions. 

The excitation temperatures derived for the lower frequency 



transition of the molecules of Table|5]lie in a range ~3-7 K, which 
is very similar to the values of T cx previously derived for CS (2- 
1) and N 2 H+ (1-0) (MGE03, MGE05). Most of the lines are also 
relatively optically thin, with a few exceptions. 

4.1.3 Results of the RADEX analysis: column density 

H2CO and SO column densities are lowest at CL1, by factors of 
~ 3.0-3.6 and ~ 1.4-2.8, respectively, compared to the other po- 
sitions, while 7V(H 13 CO + ) is the largest at CL1 by a factor of ~ 2. 
Both ICLN and ICLS show rather similar values for the column 
densities of the three molecules. 

The column density variations for CCH, C-C3H2 CN, and NO 
follow similar trends, with small differences. The CL1 position 
shows the largest column densities for all these molecules, while 
ICLS has the second largest column density. Column densities in 
CL6 tend to be slightly larger or similar to the ones in ICLN. The 
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Table 5. Physical parameters for the rest of the molecules observed in the 
four positions of LDN 673 obtained from the best-fit models after fixing 
Tk, n(H2), and 9 S to the values reported in Table[4] 



or from the hyperfine structure fitting (see Sect. IB2t . 



Position 


Species 


N(X) a 


Transition 


b 








(cm- 2 ) 




(K) 




CL1 


C-C3H2 


4.5 X 10 12 


2l, 2-1(5,1 


6.5 


0.9 




HCN 


1.5 x 10 13 


3-2 


4.9 


6.9 




S0 2 


3.7 x 10 12 


31,3-20,2 


4.5 


0.2 


CL6 


C-C3H2 


1.3 x 10 12 


2l,2-lo,l 


3.4 


0.6 


ICLN 


c-C 3 H 2 


6.0 x 10 11 


2l, 2-1(5,1 


3.3 


0.3 


ICLS 


C-C3H2 


1.9 x 10 12 


2l,2-lo,l 


3.2 


0.9 




S0 2 


5.3 x 10 12 


3l,3-2o,2 


3.2 


0.6 


Hyperfine structure calculation 


CL1 


CCH 


1.7 x 10 13 


1-0 


5.7 


0.3 




CN 


3.4 x 10 13 


1-0 


4.3 


5.9 




NO 


5.3 x 10 14 


3 5 13 
2'2 2'2 


5.4 


0.4 


CL6 


CCH 


4.7 x 10 12 


1-0 


4.5 


0.1 




CN 


9.1 x 10 12 


1-0 


3.3 


1.4 




NO 


2.2 x 10 14 


3 5 13 
2'2 2'2 


7.1 


< 0.1 


ICLN 


CCH 


5.3 x 10 12 


1-0 


4.1 


0.1 




CN 


7.3 x 10 12 


1-0 


3.1 


1.1 




NO 


1.7 x 10 14 


3 5 13 
2'2 2'2 


5.0 


0.1 c 


ICLS 


CCH 


7.0 x 10 12 


1-0 


3.9 


0.2 




CN 


1.1 x 10 13 


1-0 


3.0 


2.2 




NO 


2.5 x 10 14 


3 5 13 
2'2 2'2 


5.2 


0.1 c 



a Beam-averaged column densities. 

b Calculated for the lower frequency transition (the best determined in 
the observations). 

c Calculated after fitting the spectra assuming optically thin lines. 



column density of S0 2 , which was only detected in two positions, 
is also ~ 50% larger for ICLS than for CL1. 



4.2 Dust continuum emission 

We calculated the mass of the molecular gas traced by the dust 
from the continuum emission flux density dFrau et aT 1 120101) . Ta- 
bleU gives the H2 mass, and H2 averaged column and volume den- 
sity contained inside the FWHM contour of each of the four con- 
densations found in the 1.2-mm continuum map, assuming a dust 
temperature of 10 K. 

The calculated masses of the four condensations range from 
1.6 Mq for the SE condensation, which is also the smallest one 
(~ 0.05 pc), to 5.1 M Q for the CL1/SMM5 condensation. The H 2 
volume densities for the three centrally peaked condensations are 
rather similar, ~ 3 x 10 5 cm -3 , which is very close to the value 
derived for the CL1 position from the spectral line observations. 
The weaker and more diffuse emission of the SMM6 condensation 
is reflected in a volume density ~ 4 times lower, 7 x 10 4 cm -3 , 



Table 6. H2 average column density traced by the dust continuum emission 
on the positions of the spectral line observations for three different beam 
sizes. 



N(H 2 ) a 
r dust (10 21 cm" 2 ) 



Position 


(K) 


27" 


22" 


17" 


CL1 


10 


31.1 


34.8 


40.1 


CL6 


10 


10.8 


10.9 


10.2 


ICLN 


10 


9.0 


9.5 


11.2 


ICLS 


10 


6.0 


5.5 


4.3 



a Using the parameters given in note c of Ta- 
ble|4] 



slightly larger than the value derived for the CL6 position. Interest- 
ingly, the condensations associated with CL1/SMM5 and SMM3 
have similar extensions (~ 0.07 vs ~ 0.06 pc), volume densities, 
and, to a lesser extent, column densities. 

In order to estimate the abundances of the observed molecular 
species (see Sect. 14.3k we also calculated the H2 column density 
traced by the dust continuum emission at the four spectral posi- 
tions using three different angular resolutions: 27", 22" and 17", 
corresponding to the angular resolution of the spectral line obser- 
vations at 3.4, 2.7 and 2.1 mm, respectively (Table[6). We adopted 
a dust temperature value of 10 K for the four positions, given that 
because of the lack of any internal heating source and that the den- 
sities of the gas se em to be > 10 4 c m" 3 we do not expect higher 
dust temperatures (Goldsmith 2010). The column densities at the 
CL1 position are significantly larger than at the other three posi- 
tions, by factors of ~ 3^1 with respect to CL6 and ICLN, and of 
~5-10 with respect to ICLS. 

4.2.1 Comparison of dust and CS emission 

We convolved the dust map with the 2-point mosaic response of 
the BIMA telescope in order to properly compare the CS (2-1) and 
1.2 mm dust maps. Figure |7j shows the overlap of the CS (2-1) 
integrated emission from MGE05 with the convolved dust contin- 
uum emission at the same angular resolution, 15". In general the 
CS and 1.2 mm emissions are not correlated. For the CL1/SMM5 
condensation, the CS emission peak is slightly displaced with re- 
spect to the dust emission peak, by ~ 2". The western condensa- 
tion, SMM3, lies just at the limit of the BIMA primary beam, but 
it seems to be associated with the CS clump 12 (following MGE05 
notation). The other CS clumps found through the CLUMPFIND 
analysis done by MGE05 do not show any clear dust emission 
peaks. Some of these clumps in the northern part of the map could 
be associated with the diffuse dust emission arising from SMM6. 
More interesting is the existence of CS emission in the central parts 
of the map with very weak or no dust emission. 

To estimate the CS molecular abundance in the observed 
BIMA field (see Fig. [7J we first convolved the CS and dust con- 
tinuum maps with a beam size of 25", in order to measure the dust 
emission with a higher SNR in the regions where it is diffuse and 
weak. We then calculated the CS column density assuming LTE 
conditions, an excitation temperature of T ox =5 K, and a CS (2- 
1) optical depth line of 3, values derived by MGE05 from multi- 
transitional analysis of the CS and C 34 S molecules. The H2 col- 
umn densities were derived from the dust map adopting a temper- 
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1 9 h 20 m 55 s 50 s 
a (J2000) 



Figure 7. Greyscale of the CS (2-1) integrated emission (from MGE05) 
overlapped with the contour map of the 1 .2 mm dust emission. The angular 
resolution of both images is ~ 15". The contours of the dust continuum 
map start at 14 mJy beam - 1 with increments of 7 mJy beam -1 ; the con- 
tours of the CS (2-1) map begin at 0.496 K with increments of 0.248 K. 
The two maps were convolved with the primary beam response of the BIMA 
two-point mosaic observations presented by MGE05. The inner and outer 
dashed lines show the 0.5 and 0.25 level of the BIMA primary beam re- 
sponse. The thick crosses show the four positions observed with the IRAM 
30m telescope. 

ature Td us t=10 K and dust absorption coefficient of 0.5 cm 2 g _1 . 
Figure [8] shows the resulting CS abundance overlapped with the 
dust emission at the same 25" angular resolution. The CS abun- 
dance changes about one order of magnitude, from X(CS)~ 1CP 9 
around the CL1/SMM5 clump to A(CS)~ 1(T 8 about 1' north 
of CL1/SMM5, closer to the position of the CL6 CS clump. CS 
abundances tend to be lower at the positions of the detected dust 
condensations. 



4.3 Molecular abundances 

For each position, we calculated the molecular abundances of our 
molecules (Table IT]) from the ratio of the beam-averaged molec- 
ular column density to the corresponding beam-averaged H2 col- 
umn density from Table [6] Taking into account the uncertainties in 
the calculation of both sets of column densities, we estimate that 
the uncertainties in the determination of the molecular abundances 
should be between factors of 2-2.5. For easier comparison, we plot- 
ted in Fig. [9] the upper and lower limits that we estimated for the 
abundances of each molecular species in CL6, ICLN, and ICLS 
relative to the abundance in CL1. 

The abundances of CCH, CN, and 0C3H2 at CL1 are similar 
to the values of the other positions, but for the rest of the molecules 
the abundances are clearly lower, by a factor of ~ 2 for H 13 CO + 




1 9 h 20 m 55?0 50 s 
a (J2000) 



Figure 8. Colour image of the map of the CS abundance, overlapped with 
the contour map of the 1.2 mm continuum emission. The angular reso- 
lution of both images is 25". The two maps were convolved with the 
primary beam response of the BIMA two-point mosaic observations pre- 
sented by MGE05. The inner and outer dashed lines show the 0.5 and 0.25 
level of the BIMA primary beam response. The contours are 10%, 15%, 
25%, 35%, 55%, 75%, and 95% of the 1.2 mm continuum peak emission, 
(66.2 mJy beam -1 ). The grey crosses show the four positions observed 
with the IRAM 30m telescope. 



and NO, and of ~ 10 for SO, S0 2 , and H 2 CO. On the other hand, 
the abundances are systematically larger for all species at ICLS, but 
only by a factor of < 3 with respect to the abundances at CL6 and 
ICLN. 



5 DISCUSSION 

5.1 Fragmentation in the core 

The most intense emission in the central region of the 1.2-mm 
dust continuum map is contai ned inside the HPBW of the low- 
angular resolution NH3 map of lSepulveda et alj J201 lh . The BIMA 
observations of MGE03 also only detected N2H" 1 " emission in a 
region almost coinciding with CL1/SMM5 (the other dust contin- 
uum sources lie outside the BIMA primary beam). The dust and 
N2H + emissions in CL1/SMM5 are similarly elongated in an ap- 
proximately N-S direction with a close proximity between both 
emission peaks. A second clump, located just North of CL1/SMM5 
could also be present in the dust emission, but it is not so clearly 
seen as in the N 2 H + map. The sizes of the CL1/SMM5 conden- 
sation in the continuum and molecular (CS and N2H + ) maps are 
very similar, ~0.07-0.09 pc, but the gas mass traced by the dust 
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Table 7. Molecular abundances calculated from the ratio of the beam- 
averaged molecular column density to the beam-averaged H2 column den- 
sity derived from the dust continuum emission. 



Species 



Position Beam a 



[X]/[H 2 ] 



CCH 


CL1 


27" 


5.1 


X 


10- 10 




CL6 




1.3 


X 


10 -io 




ICLN 




5.8 


X 


10 -io 




ICLS 




1.2 


X 


10- 9 


CN 


CL1 


22" 


9.9 


X 


10- 10 




CL6 




8.3 


X 


10- 10 




ICLN 




7.6 


X 


10 -io 




ICLS 




2.0 


X 


10- 9 


C-C3H2 


CL1 


27" 


1.4 


X 


10 -io 




CL6 




1.2 


X 


10- 10 




ICLN 




6.7 


X 


10- 11 




ICLS 




3.2 


X 


10 -io 


HCN 


CL1 


12" 


3.0 


X 


10 -io 


H2CO 


CL1 


17" 


2.1 


x 


10 -io 




CL6 




2.7 


X 


10- 9 




ICLN 




2.7 


X 


10- 9 




ICLS 




5.8 


X 


10- 9 


Til 3/-v-\-t- 


CL1 


(1-7// 

27 


2.7 


X 


1 n — 11 

10 




CL6 




1.2 


X 


10- 11 




ICLN 




4.4 


X 


10- 11 




ICLS 




6.6 


X 


10- 11 


NO 


CL1 


17" 


1.3 


X 


10-s 




CL6 




2.2 


X 


10-8 




ICLN 




1.5 


X 


ID" 8 




ICLS 




5.7 


X 


10-8 


SO 


CL1 


22" 


1.4 


X 


10- 9 




CL6 




1.2 


X 


10-8 




ICLN 




6.9 


X 


10- 9 




ICLS 




1.1 


X 


ID" 8 


so 2 


CL1 


22" 


1.1 


X 


10 -io 




ICLS 




9.6 


X 


10 -io 



a Beam-size over which the molecular and H2 col- 
umn densities were averaged. 



is ~ 4 times larger than the masses estimated, assuming standard 
abundances, from the CS and N2H + maps by MGE05. 

We showed in Sect. I4.2.T1 that the CS and dust emissions do 
not seem to be particularly correlated (see Figs.|7]and|8]: the most 
intense CS emission around CL1/SMM5 follows the edges of the 
most intense dust emission and many of the clumps detected in the 
combined FCRAO and BIMA maps are located in regions with no 
or very weak dust emission. These results suggest that there is some 
low density/small-size structure which is revealed by the CS lines, 
but it is not detected at 1.2-mm, probably because the sensitivity of 
our observations is not high enough or due to the spatial filtering of 
the chopping bolometer observations. 

The SMM8 dust condensation is also coincident with rela- 
tively weaker H 13 CO+ and CS emission detected with the FCRAO 
telescope by MGE05 at the tip of a very weak NW-SE filamen- 
tary structure or thin ridge found in the molecular line maps . This 
filamentary structure is undetected in the dust continuum map. 

In order to study the fragmentation in the core, we calculated 
the Jeans-length corresponding to the densities of the two main 
dust condensations, CL1/SMM5 and SMM3, following Hartmann 




CCH CN c-C 3 H 2 H 2 CO H 13 CO + NO SO S0 2 
molecular species 

Figure 9. Relative abundances of the molecules detected at the CL6 (cir- 
cles), ICLN (inverted triangles), and ICLS f squares) positions with respect 
to the abundance at CL1. The error bars indicate the estimated uncertainties 
in the determination of the molecular abundances. 



dl998h : 

Aj = 



na 
~~Gp 



1/2 



ai9pc (i^) 



1/2 



riH, 



10 4 cm- 



-1/2 



(1) 



with the sound speed a, the gravitational constant G, and the den- 
sity p. The separation between the condensations is ~ 50", which 
at the distance of LDN 673 corresponds to ~ 0.073 pc (~ 15000 
AU). We used the average of the densities estimated for both con- 
densations from the continuum observations, ~ 3.3 X 10 5 cm -3 , 
as the volume density, riH 2 . and assumed T = 10 K. We obtained 
a Jeans-length value of ~ 0.033 pc, about 2 times smaller than 
th e measured separatio n on the sky, and similar to the values found 
bv lTsitali et alj OOlOh for some other condensations in the region. 
This value of the Jeans-length is also smaller than the sizes of both 
condensations, which is an indication that they may still be sub- 
ject to more thermal fragmentation. Indeed, there is evidence of 
fragmentation in the CL1/SMM5 condensation, since this clump 
is associated with two N2JT 4 " condensations separated by ~ 30" 
(MGE05), which corresponds to ~ 0.044 pc (~9000 AU), very 
close to the Jeans-length previously estimated. On the other hand, 
most of the CS clumps detected by MGE05 and not traced by the 
1.2-mm dust map have typical sizes of 0.04 pc and densities of a 
few times 10 4 cm -3 . For a density of 7 x 10 4 cm -3 (the CL6 
density), the Jeans-length is 0.07 pc. Thus, these clumps are un- 
likely to have been formed as a result of thermal fragmentation and 
we would need to invoke non-thermal pro cesses (turbulence, MHD 
waves, etc) to explain their existence 1 Falle & Hartquist 20021; 
Vazquez-Semadeni et al 1 l2005l ; I Klessen et a l. 2005; V an Loo et al .1 
200 8|). This would agree with the transient nature of this clumps 
postulated by MGE05. 



5.2 The nature of the clumps and the inter-clump medium 

The interpretation of single-point observations of molecular lines as 
the ones presented here have to be taken with caution, because they 
will never provide the same level of detail of the chemistry as a full 
map would. However, our choice of positions with presumed dif- 
ferent physical conditions, coupled with the dust continuum map, 
provides us with a sample of molecular abundances that can show 
some of the properties of our core. 
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MGE05 already showed that CS was probably depleted in 
CL1, while N2H + was clearly centrally condensed. The estima- 
tion of the CS abundance in Fig. [8] confirms this result for the 
SMM5, SMM3, and SMM6 dust condensations. Additionally, Ta- 
ble|7]shows that H 2 CO, SO, SO2, and, to a lesser extent, H 13 CO+ 
are also depleted for CL1 with respect to the abundances at the 
other positions. On the other hand, the abundances of CCH, c- 
C3H2, and CN in CL1 do not seem to differ much from the values in 
CL6 and ICLN. Finally, the NO abundance at CL 1 is also ~ 40% 
lower than the one at CL6, and similarly to what lAkvilmaz et al.l 
d2007h found in LI 544, the NO molecule could be partially de- 
pleted in CL1, while N2H + is not; 

The or i ginal mode ls of iTavlor et al.l d 19961 1 19981) and 
iGarrod et al l d2005l 120061) required the action of depletion to ob- 
tain the distribution of molecular abundances that would explain 
the results of MGE03 and MGE05. As has been extensively found 
in different starless cores, 'early-time', carbon-bearing molecules, 
would be depleted in the more central regions of evolved starless 
cores, wherein 'late-time', nitr ogen-bearing molecules, would still 
be present in the gas (see e. g JCaselli et al.ll999l : lHilv-Blant et all 
l20ld : lTafalla et alj|2006h . Thus, the results for CS, N 2 H+, H 2 CO, 
SO, and H 13 CO + are approximately as expected, given the density 
of the CL1/SMM5 condensation. However, the low or n o depletion 
of CC H, C-C3H2, and CN is much more unexpected. iFrau et all 
j2012l) also found a similar result in several cores of the Pipe Neb- 
ula, where SO shows depletion, but 0C3H2 and CCH are more 
intense in the cores with higher Av. 

CL6 and the inter-clump positions show a poorer variety of 
molecules compared to that of CL1. Additionally, there is little evi- 
dence of depletion of the detected molecules - the abundances of 
early- time molecules are usually the largest in the inter-clump po- 
sitions - with no detection of N 2 H + . All this seems to indicate that 
the gas in CL6, ICLN, and ICLS is chemically young. The gas den- 
sity in CL6 is more similar to the densities found in the inter-clump 
medium (see Sect. 14. L2) than to the density in CL1. However, the 
number density derived for CL6 is still a factor of ~ 3.5 larger 
with respect to the two inter-clump positions, which indicates that 
the physical conditions are probably different. This factor is signifi- 
cant and probably a lower limit to the true density contrast, because 
the probable contribution to the emission from nearby clumps will 
have the effect of measuring a density higher than the true inter- 
clump density in ICLN and ICLS (see Sect |4.1.lT ). In any case, it 
seems clear that the inter-clump positions are denser than initially 
expected, probably from a few times 10 3 to 10 4 cm~ 3 . 

5.2.1 Comparison with models of transient clumps 

In order to have a more quantitative comparison, we compared the 
results of our observations to the "standard" model of IGarrod et al .1 
j2005ll2006h . We found that this model is able to explain reason- 
ably the chemistry of the CL6 and ICLN positions. However, we 
used a "high-density" model th at reached a higher peak density, 
~ 5.5 x 10 5 cm" 3 , in 10 6 yr dGarrod 120051) to match better the 
chemistry in CL1. We did not compare the models to ICLS due to 
the contamination issues discussed in Sect. 14. 1.71 We consider that 
the real evolutionary properties of the gas we observed are pro- 
bably between these two models and some scaling of abundances 
would be needed to really fit our results, which is outside the pur- 
pose of this discussion. Given the traditional difficulty of matching 
the values of chemical models to observations, we looked for the 
time ranges in which the abundances of most of the molecules in 
the models (> 70%) were within a factor of 3 from the observed 



values, in order to account for the multiple sources of uncertainties 
in both observations and models. 

The molecular abundances in CL1 are compatible with a 
clump in the "high-density" model in a time range of 3-7 x 10 J yr. 
HCO + is the main outlier, depending on the 12 C/ 13 C abundance 
ratio used, and "late-time" molecules, NO and SO, tend to push 
the time range to later times. The abundances in CL6 are better 
reproduced with a clump following the standard model at times ei- 
ther 5-9xl0 5 yr or 1.4-1.6 x 10 6 yr. So, the clump in CL6 could 
be very close to the peak density maximum or just past it in the 
dissipation phase. We also found that the results for ICLN could 
be compatible with a clump of 0.9-1. 6 xlO 6 yr for the standard 
model. This could explain the relatively low contrast in density 
between CL6 and ICLN, if we think that the gas in ICLN could 
be the remains of a clump that has almost dispersed and what we 
found was lower density gas enriched by the molecular abundances 
of previous s t ages. This would reinforce one of the predictions of 
IGarrod et al.l d2006l) models: that the formation and dispersion of 
transient clumps would result in a chemistry of the gas generally 
"young", except in the denser and more evolved clumps, and a gas 
chemical composition that would show signs of "enrichment" of 
abundances of some molecular species in the lower density mate- 
rial. 

IGarrod et al.ld2006l) also showed that a distribution of transient 
clumps at different stages of evolution (different points in time) 
could qualitatively reproduce the medium- and high-angular reso- 
lution maps of MGE03 and MGE05. Our results here suggest that 
those clumps could not only be in different stages of evolution but 
also reach different peak densities, i.e. follow different evolutionary 
patterns. The mix of clumps with different evolutionary properties 
can shed a different light on the discussion about the molecular de- 
pletion between different clumps. We found that, for the times best 
fitted by the models, the differences in the abundances of CCH, c- 
C3H2, and CN between the "high-density" and "standard" models 
are considerably smaller than for later-time molecules, such as NO. 
This would simply explain the apparent low or no depletion of the 
former molecules compared to the latter. 



6 CONCLUSIONS 

We studied the properties of the clump and inter-clump gas in 
the CS starless core of LDN 673 using the 30-m IRAM telescope 
through the emission of several spectral lines in the 3-, 2-, and 1- 
mm bands, in four positions associated with already detected, but 
physically and chemically different, clump gas (positions CL1 and 
CL6), or with inter-clump gas (positions ICLN and ICLS). We de- 
tected 19 spectral transitions of 10 molecular species at least in 
one of the four positions. We complemented the spectral observa- 
tions with the mapping of the 1.2-mm dust continuum emission in 
LDN 673, which allowed us to obtain a more reliable estimate of 
the volume density of the gas and of the abundances of the detected 
molecules. The main results of our study are: 

(i) The dust continuum observations revealed four emission con- 
densatio ns in the region, ro ughly coinciding with SMM sources 
found bv lVisser et al.l yOOZ), three of them with a round shape and 
centrally peaked emission: CL1/SMM5 is the most intense and co- 
incides with the N 2 H+ clump found by MGE03; SMM3, located 
~ 50" to the West, was not previously found in the molecular 
observations; and SMM8, located ~ 2'5 SE of the centre of the 
map, and coinciding with an emission enhancement in the C 34 S and 
H 13 CO+ maps of MGE05. Finally, there is a more diffuse emission 
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strip extended E-W ~ 2' north of CL1/SMM5. The masses of these 
condensations range from 1.6 to 5.1 Mq and sizes 0.05-0.07 pc. 
The northern condensation is sensibly more diffuse. 

(ii) We made a radiative transfer analysis using the RADEX 
code to determine the molecular column densities, volume density, 
and kinetical temperature at each position. The best fits to our ob- 
servations found T K ~ 10 K for the CL1 and CL6, and ~ 20 K for 
the ICLN and ICLS positions; and volume densities ranging from 
~ 3.6 x 10 5 cm" 3 at CL1 to ~ 1.7 x 10 4 cm" 3 at ICLS. 

(iii) CL1 presents the largest column densities for most of the 
molecules (H 13 CO+, c-C 3 H 2 , CCH, CN, and NO), but the smaller 
column densities for H2CO and SO. The estimated molecular abun- 
dances of most of the molecules, are smallest at CL1, while the 
abundances at CL6 and ICLN tend to be very similar. 

(iv) The comparison of the dust continuum emission to previous 
interferometer and single-dish observations shows very little or no 
correlation between the CS (2-1) and the 1.2-mm emission, with 
most of the clumps found in the CS emission not detected in the 
continuum observations. On the other hand, the dust continuum and 
N2H + emissions are much more similar. 

(v) We found that the condensations detected in the dust con- 
tinuum map may be subject to more thermal fragmentation, which 
might be already happening in CL1/SMM5. However, most of the 
clumps detected in CS, but undetected in the dust continuum ob- 
servations, are unlikely to have been formed by thermal fragmenta- 
tion, and would need other mechanisms to explain their formation. 
These clumps were proposed to be transient by MGE05. 

(vi) The chemistry of CL1 appears to be much more evolved 
than for the other positions, with signs of depletion of several 
molecular species (CS, H 2 CO, SO, H 13 CO+, and partially NO), 
but also relatively unexpected high abundances of CCH, 0C3H2, 
and CN. The density contrast between CL6 and the two inter-clump 
positions, which are denser that initially expected, is relatively low. 
The gas at the CL6 and inter-clump positions seems to be generally 
chemically "young". 

In summary, the central condensation, CL1/SM M5, is proba- 
bly approach ing the 'peak- time' state of the models of lGarrod et all 
J2005l . l2006l) and it is also the place with a larger probability of fu- 
ture undergoing of star formation (see MGE05). The SMM3 con- 
densation is probably in a very similar state, from its shape and gas 
density, but we cannot firmly determine it until further observations 
provide us with more information about its chemistry. The low den- 
sity contrast between CL6 and the inter-clump positions and their 
similar young chemical age seems to support the idea of the pres- 
ence of lower dens ity transient clumps in the core as proposed by 
lGarrodetal](l2006l) . 
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APPENDIX A: 
LINES 



LINE PARAMETERS OF THE OBSERVED 



Tables [ATI to I A4I show the line parameters obtained using a Gaus- 
sian fit for the transitions detected in the CL1, CL6, ICLN, and 
ICLS positions, respectively. Each table lists the molecular species 
and transition and the four parameters resulting from the Gaussian 
fit: line intensity (in main-beam units), central velocity of the line, 
line-width, and intensity integrated under the Gaussian. 

Table IA5I shows the hyperfine structure fit parameters, ob- 
tained using the HFS method of the CLASS package, for the CN, 
CCH and NO transitions detected at 3-mm. The table lists: posi- 
tion, molecular species and transition, Ar, central velocity of the 
hyperfine line of reference, and optical depth of this line. 

Finally, Table I A6l gives the 3cr upper limits for the line inten- 
sities of the transitions not detected at each position. 



APPENDIX B: DETERMINATION OF THE BEST-FIT 
SOLUTIONS USING THE RADEX MODELLING 

Bl Method for multiple-line detected molecules 

We used for our calculations the line intensities, T^, of H2CO, 
H 13 CO + , and SO, which have at least three observed transitions at 
each of the four positions. We ran a grid of RADEX models using 
the following ranges in temperature, volume density, and molecu- 
lar column density: T K = 10-30 K; n(H 2 ) = 10 3 -10 7 cirT 3 ; and 
N(X) = 10 11 -10 15 cm -2 . For each set of temperature, volume 
density, and column density values, RADEX provides the line in- 
tensity for each i-transition, T,l x . The expected, or calculated, mean 
beam temperature is given by 



T' 



Tii i 
rxVbf 



(Bl) 



where rjl * is the beam filling-factor of the i-transition. Assuming a 
2-D Gaussian distribution for the source emission profile, the filling 
factor can be expressed as 



We obtained the "best-fit mo del" after searching for the min- 
imum of the reduced \ 2 function jNummelin et aljfeoOOl) . re- 
sulting from comparing the measured and calculated line intensities 
and line intensity ratios for the three aforementioned molecules: 



xl 



E 



(B3) 



R r 



T 1 

J-ni 



T 

1 r, 



where T^ s is the observed line intensity of the lowest frequency 
transition (which is usually the best determined) of the k molecule, 
R^ hs is the observed line intensity ratio between the i- and j transi- 
tions of the k molecule, o]£ s is the Ict uncertainty associated to the 
observed line intensity of the lowest frequency of the k molecule, 
a oli k ' s m e uncertainty associated to R% bs , n is the number of data 
points and p is the number of free parameters. The calc sub-indexes 
indicate the equivalent values for the line intensities derived with 
RADEX. For the undetected transitions, we only included in our 
calculations the intensities from RADEX that were larger or equal, 
within errors, than the observed ones. 

Since the real size of the emitting area at each position was 
unknown, we used two sets of values of filling-factors in our analy- 
sis: r)lf ~ 1, obtained by choosing 9 S — 1000", and r\ b j < 1. For 
the second case, we tested values between 30", similar to the beam 
size at 3-mm, and 50", the deconvolved size of the CL1/SMM5 
dust condensation as well as the size of the CL6 clump derived by 
MGE05. 



B2 Method for the rest of molecules 

Following a method similar to the one described in Sect. IB II but 
fixing the values of n(H.2), 7k, and 8 S to the best-fit values de- 
termined for each position, we looked for the best set of values of 
column density, N, excitation temperature, T ex , and line opacity, r 
for the rest of molecules. 

We did not follow this method for the CCH lines, because 
the collisional rates for this molecule are still unknown (see e.g., 
IPadovani et~al]|2009l) . nor for the CN and NO transitions, because 
the RADEX code does not take into account the hyperfine line 
structure, which in this case ends up underestimating the opacity 
of the lines. Instead, we calculated T cx and r of the lines from the 
simultaneous fitting of the hyperfine components. In this last case, 
we assumed abeam-filling factor, r\ = 1. 



Vbf 



+ 



Oi 2 ' 
rnb 



(B2) 



where 6 % mb is the FWHM of the antenna main beam at the frequency 
of the i-transition, and 9 S is the FWHM source size. 
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Table Al. Line parameters for the transitions detected in Clump 1 (CL1) 



Tub Vlsr Line area 



Transition 




(K) 


(kms -1 ) 


(kms- 1 ) 


(Kkms" 


- 1 ) 


CCH 2-1 





35 ± 0.07 


±7.10 ±0.04 


0.48 ± 0.09 


0.18 ±0 


03 


CN2-1 





79 ± 0.08 


±7.00 ±0.01 


0.41 ± 0.03 


0.35 ±0 


02 


C-C3H2 2 ,l-ll,2 


1 


46 ± 0.02 


±7.16 ±0.01 


0.49 ± 0.01 


0.76 ±0 


01 


C-C3H2 4 M -3 ,3 





52 ± 0.06 


±7.15 ±0.01 


0.51 ± 0.03 


0.28 ±0 


02 


HCN 3-2 





88 ± 0.26 


±6.84 ±0.03 


0.28 ± 0.08 


0.27 ±0 


06 


H 2 CO 2l >2 -ll,l 


1. 


73 ± 0.06 


±7.02 ±0.01 


0.54 ±0.02 


1.00 ±0 


02 


H 2 C0 3i, 3 -2i,2 


1 


14 ±0.07 


±7.00 ±0.01 


0.46 ± 0.02 


0.56 ±0 


03 


H 2 C0 3i, 2 -2i,i 





75 ± 0.06 


±7.03 ±0.01 


0.41 ± 0.03 


0.33 ±0 


02 


H 13 CO+ 1-0 


1. 


44 ± 0.04 


±7.23 ±0.01 


0.52 ± 0.01 


0.79 ±0 


02 


H 13 CO+ 2-1 


1. 


04 ±0.15 


±7.19 ±0.02 


0.38 ± 0.06 


0.42 ±0 


.05 


H 13 CO+ 3-2 





39 ± 0.13 


±7.25 ±0.03 


0.23 ± 0.09 


0.10 ±0 


03 


NO 5 1-3 6 

w 2'2 2'2 





39 ±0.12 


±7.14 ±0.06 


0.50 ±0.11 


0.21 ±0 


.05 


SO 3 2 -2i 


3 


45 ± 0.03 


±7.12 ±0.01 


0.52 ± 0.01 


1.89 ±0 


01 


SO 4 5 -3 4 





17 ± 0.04 


±7.07 ±0.03 


0.27 ±0.06 


0.05 ±0 


01 


SO 65-54 





75 ± 0.09 


±7.03 ±0.01 


0.35 ± 0.02 


0.28 ±0 


02 


SO2 31,3-20,2 





19 ± 0.03 


±7.03 ±0.03 


0.42 ± 0.06 


0.08 ±0 


01 



Table A2. Line parameters for the transitions detected in Clump 6 (CL6) 







Tmb 


Vlsr 


AV 


Line area 


Transition 




(K) 


(kms" 1 ) 


(kms" 1 ) 


(Kkrns" 1 ) 


C-C3H2 2 , 


l-ll,2 


0.20 ±0.01 


±7.23 ±0.01 


0.46 ± 0.03 


0.10 ±0.01 


H 2 CO 2l,2 


-ll.l 


1.25 ± 0.07 


±7.16 ±0.01 


0.67 ±0.02 


0.90 ±0.03 


H 2 C0 3i, 3 


-2l,2 


0.77 ±0.06 


±7.20 ±0.01 


0.56 ± 0.04 


0.45 ±0.02 


H 2 C0 3i, 2 


-2i,i 


0.38 ± 0.05 


±7.17 ±0.02 


0.48 ± 0.06 


0.19 ±0.02 


H 13 CO+ 1- 


-0 


0.51 ± 0.03 


±7.35 ±0.01 


0.60 ±0.03 


0.32 ±0.01 


H 13 CO+ 3 


-2 


0.35 ± 0.10 


±7.35 ±0.43 


0.26 ± 0.12 


0.10 ±0.04 


SO 3 2 -2i 




2.87 ±0.04 


±7.24 ±0.01 


0.52 ± 0.01 


1.60 ±0.02 



Table A3. Line parameters for the transitions detected in Inter-clump North (ICLN) 







Tmb 


Vlsr 


AV 


Line area 


Transition 




(K) 


(kms" 1 ) 


(kms- 1 ) 


(Kkrns- 1 ) 


c-C 3 H 2 2 , 


l-ll,2 


0.10 ±0.02 


±7.27 ±0.02 


0.40 ±0.05 


0.04 ±0.01 


H 2 CO 2i, 2 - 


-ll.l 


0.79 ± 0.07 


±7.19 ±0.02 


0.95 ± 0.04 


0.80 ±0.03 


H 2 C0 3i, 3 


-2l,2 


0.37 ±0.06 


±7.32 ±0.03 


0.79 ± 0.08 


0.31 ±0.03 


H 2 C0 3i, 2 


-2i,i 


0.14 ±0.04 


±7.34 ±0.06 


0.84 ±0.14 


0.12 ±0.02 


H 13 CO+ 1- 


-0 


0.35 ± 0.02 


±7.34 ±0.01 


0.61 ± 0.03 


0.23 ±0.01 


SO 3 2 -2i 




1.74 ±0.03 


±7.23 ±0.01 


0.71 ± 0.01 


1.31 ±0.01 


SO 6 5 -5 4 




0.18 ±0.06 


±7.51 ±0.07 


0.50 ±0.17 


0.10 ±0.03 
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Table A4. Line parameters for the transitions detected in Inter-clump South (ICLS) 





Tmb 


^LSR 


AV 


Line area 


Tr n n q i ti n n 

1 1 Ull>llHMI 




(km s — 1 ) 


(km s —1 ) 


IK km i — ^ 


CN 1-0 


0.26 ±0.03 


±7.20 ±0.10 


0.65 ± 0.10 


0.10 ±0.01 


C-C3H2 2 ,i-li, a 


0.25 ± 0.03 


±7.27 ±0.02 


0.37 ±0.04 


0.10 ±0.01 


H 2 CO 2l, 2 -ll,l 


0.90 ±0.07 


±7.21 ±0.01 


0.46 ± 0.03 


0.44 ± 0.02 


H 2 CO 3i,3-2i, 2 


0.37 ±0.05 


±7.23 ±0.02 


0.64 ± 0.06 


0.25 ±0.02 


H 2 CO 3i, 2 -2i,i 


0.33 ± 0.08 


±7.24 ±0.02 


0.38 ± 0.06 


0.14 ±0.02 


H 13 CO+ 1-0 


0.38 ± 0.02 


±7.40 ±0.01 


0.44 ± 0.02 


0.18 ±0.01 


H 13 CO+ 2-1 


0.24 ±0.07 


±7.44 ±0.06 


0.36 ± 0.15 


0.09 ± 0.03 


SO 3 2 -2i 


1.99 ± 0.03 


±7.32 ±0.01 


0.47 ±0.01 


1.00 ±0.01 


SO2 3i,3-2q,2 


0.14 ±0.04 


±7.23 ±0.03 


0.23 ± 0.07 


0.04 ±0.01 



Table A5. Line parameters for the transitions with hyperfine structure 











At 




AV 




Position 


Transition 




(K) 


(km s" 1 ) 


(km s" 1 ) 


r 


CL1 


CCH 1-0 






0.90 ±0.07 


±7.19 ±0.01 


0.45 ± 0.02 


0.33 ±0.29 




CCH 2-1 






0.88 ±0.72 


±7.13 ±0.02 


0.28 ± 0.09 


1.97 ± 2.75 




CN 1-0 






7.76 ±0.83 


±7.09 ±0.01 


0.35 ± 0.01 


5.90 ±0.83 




CN2-1 






0.83 ±0.18 


±7.00 ±0.01 


0.40 ±0.03 


0.13 ±2.31 




NO §,§- 


1 
2 


:i 
2 


0.78 ±0.11 


±7.09 ±0.01 


0.48 ± 0.04 


0.38 ±0.28 


CL6 


CCH 1-0 






0.16 ±0.01 


±7.24 ±0.02 


0.58 ± 0.05 


< 0.10 




CN 1-0 






0.58 ±0.07 


±7.24 ±0.01 


0.57 ±0.04 


1.37 ±0.38 




NO f ,§- 


1 
2 


3 
2 


0.36 ± 0.03 


±7.21 ±0.02 


0.54 ±0.07 


< 0.10 


ICLN 


CCH 1-0 






0.12 ±0.01 


±7.29 ±0.02 


0.76 ± 0.07 


< 0.10 




CN 1-0 






0.28 ±0.05 


±7.22 ±0.02 


0.62 ± 0.05 


1.09 ±0.47 




NO §,§- 


1 
2 


3 
2 


0.17 ±0.02 


±7.20 ±0.04 


0.64 ±0.10 


0.10 a 


ICLS 


CCH 1-0 






0.22 ±0.06 


±7.32 ±0.02 


0.50 ±0.05 


0.22 ±0.65 




NO §,§- 


1 
2 


:i 
2 


0.18 ±0.03 


±7.26 ±0.07 


0.89 ± 0.22 


0.10 a 



a Obtained after assuming optically thin lines. 
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Table A6. Line intensity 3cr upper limits for the transitions not detected 



Position 

CL1 CL6 ICLN ICLS 

Tmb Tmb Tmb Tmb 



Transition 


(K) 




(K) 




(K) 


(K) 


CCH 2-1 




< 


0.15 


< 


0.14 


< 0.20 


CCH 3-2 


< 0.72 


< 


0.58 


< 


0.50 


< 0.72 


CN2-1 




< 


0.14 


< 


0.12 


< 0.12 


CS 5^1 


< 0.19 


< 


0.34 


< 


0.28 


< 0.29 


C-C3H2 41,4-30,3 




< 


0.13 


< 


0.17 


< 0.19 


HCN 3-2 




< 


0.60 


< 


0.52 


< 0.49 


H 13 CO+ 2-1 




< 


0.24 


< 


0.20 




H 13 CO+ 3-2 








< 


0.36 


< 0.37 


NO - --- - 

2'2 2'2 




< 


0.43 


< 


0.31 


< 0.59 


SO 4 5 -3 4 




< 


0.10 


< 


0.10 


< 0.22 


SO 6 5 -5 4 




< 


0.31 






< 0.24 


SO 7 6 -6 5 


< 0.87 


< 


0.45 


< 


0.42 


< 0.49 


SO2 31,3-20,2 




< 


0.14 


< 


0.12 




SO2 7i, 7-60,6 


< 0.26 


< 


0.35 


< 


0.25 


< 0.30 
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